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Determination of the Coronal Magnetic Field by Hot Loop 
Oscillations Observed by SUMER and SXT 

Tongjiang Wang 1 , Davina E. Innes 2 , and Jiong Qiu 1 
ABSTRACT 

We apply a new method to determine the magnetic field in coronal loops 
using observations of coronal loop oscillations. We analyze seven Doppler shift 
oscillation events detected by SUMER in the hot flare line Fe XIX to obtain 
oscillation periods of these events. The geometry, temperature, and electron 
density of the oscillating loops are measured from coordinated multi-channel soft 
X-ray imaging observations from SXT. All the oscillations are consistent with 
standing slow waves in their fundamental mode. The parameters are used to 
calculate the magnetic field of coronal loops based on MHD wave theory. For the 
seven events, the plasma f3 is in the range 0.15—0.91 with a mean of 0.33±0.26, 
and the estimated magnetic field varies between 21—61 G with a mean of 34±14 
G. With background emission subtracted, the estimated magnetic field is reduced 
by 9%— 35%. The maximum backgroud subtraction gives a mean of 22±13 G 
in the range 12—51 G. We discuss measurement uncertainties and the prospect 
of determining coronal loop magnetic fields from future observations of coronal 
loops and Doppler shift oscillations. 

Subject headings: Sun: corona — Sun: flares — Sun: oscillations — Sun: UV 
radiation — Sun: X-rays, gamma rays 



1. Introduction 

Magnetic field governs the structure and dynamics of the Sun's corona. Although the 
structure and topology of the coronal magnetic field can be indirectly traced in imaging 
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observations of coronal loops in X-rays and EUVs, a direct measurement of the coronal 
magnetic field remains a very difficult and challenging problem. In the past, efforts have 
been made with microwave, optical and EUV observations to measure coronal magnetic fields. 
The microwave gyroresonance magnetometry technique is used to measure strong (more than 



a few hundred gauss) active region field str engths in the low corona (e.g. iGary fe Hurford 



1994 ; ISchmelz et al. 111994 ; iLee et al. Ill999l ). and Faraday rotation observations of polarized 
radiation by radio sources can be used to measure the mean line-of-sight coronal field (e.g. 
Sakurai fe Spangler Ill994l ). At optical and EUV wavelengths, the Zeeman and Hanle effects 
are used to determine the coro nal magnet i c field, but their applications are limited because 



the signals are weak. Recently, iLin et al. I (120001 . 120041 ) accurately measured the line-of-sight 
component of the coronal magnetic field using Zeeman splitting observations in the near- 
infrared line Fe XIII 10747 A and obtained a field strength of 4—33 G at a height ranging 



betwe en 0.10 and 0.15 of a solar radius off the limb above active regions. iRaouafi et al 



(120021 ) made the first attempt to measure the coronal field from the linear polarization of 
the O VI 1032 A line due to the Hanle effect and deduced a coronal field of about 3 G in a 
polar hole. 



A new m ethodology, known as MHD coronal seismology (e.g. iRoberts et al. 1 11984 



Roberts 1120041 ). has bee n developed to diagnose plasm a properties of coronal loops exhibit- 
ing oscillation patterns. iNakariakov fc Ofman I (120011 ) first applied this method to TRACE 
observations of global fast kink-mode loop oscillations in two events, yielding a loop field 
strength in the range of 4—30 G. Taking into consideration the plasma d ensity in one event. 



they r efined their determination to 13±9 G. With the same method, lAschwanden et al. 



( 120021 ) obtained the magnetic field strength ranging between 3 and 90 G for 26 oscillating 
loops observed by TRACE. 

In this paper we analyze slow mode MHD waves in coronal loops to determine their 
magnetic field. Standing slow mode waves, triggered by small flares, were recently discovered 
in spectroscopic observations taken by SUMER on- board SOHO, by t heir D oppler shift 



opp 

signatures in the flare lines Fe XIX and Fe XXI (e.g IWang et al.l 120021 . l2003al ibl). Similar 



strongly damped Doppler shift oscillations were also detected by BCS on Yohkoh m many 



flares , whic h were also interpreted as signatures of standing slow mode waves (IMariska 



20051 . 120061 ). The plasma (3 in hot (> 6 MK) loops can be much larger than the expected 
mean value (~0.01) in the corona, which implies a relatively strong coupling between the 
magnetic and thermal perturbations. Therefore, the loop field strength may also be deduce d 
from diagnostics of standing slow mode oscillations in hot coronal loops (IRoberts I I2004J ) . 
In this study we will test this idea. In the following text, the observations are described in 
Sect. 2. The method for deducing the magnetic field strength is given in Sect. 3. We present 
the results and uncertainty analysis in Sect. 4, and the conclusions and discussions in Sect. 5. 
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Observations 



Spectral observations of active region coronae were obtained by SUMER in sit-and-stare 
mode. Two neighboring active regions on the east limb were observed with a cadence of 90 s 
using the 300" x 4" slit during 2000 September 16-20. Five spectral lines, Fe XIX A 1118 (6.3 
MK), Ca XV A 1098 and A 555x2 (3.5 MK), Ca XIII A 1134 (2.2 MK), and Si III A 1113 (0.06 
MK) were transmitted, with a 2.2A wide window for each line. After processing the raw 
data following standard procedures (decompression and corrections of flat-field, detector 
distortions, deadtime, and gain effects), a single Gaussian was fit to each line profile to 
obtain a Doppler shift time series at each spatial pixel (e.g. Fig. [T]c). A large number of 
flarelike brightenings were revealed in the hot flare line, Fe XIX. They have an occurrence 
rate 3—14 per hour and their lifetimes range from 5—150 min, with an average of about 25 
min (jWang et al. II2006I ). From about 180 events identified by the visual inspection, more 
than 20% of them are found to be associated with Doppler shift oscillations. 

For six Doppler shift events in this dataset, coordinated imaging observations were ob- 
tained by SXT on Yohkoh. In addition, coordinated SXT observations were also made for 
an oscillation event reported in an earlier study by IWang et al.l (120021 ). thus this event is 
also included in this study. Analyses of joint SXT and SUMER observations of coronal loop 
oscillations allow us to identify the oscillating loops, determine their geometric parameters, 
and derive their temperature and density. In order to compare SXT and SUMER observa- 
tions, the scale of the SXT full (2". 5 pixel -1 ) and half-resolution images is increased by a 
factor of about 1% to account for the different orbits of Yohkoh and SOHO. Then the SXT 
images are coaligned with the SUMER slit according to their pointing information. The 
Fe XIX brightenings along the slit coincide accurately with the slit intersection on the SXT 
loops, indicative of a good coalignment (see Figs. [B15]a and [I]|5]&) . The average temperature 
and emission measure of the loops are derived from the SXT images using the filter ra- 
tio met hod implemented in the standard SolarS oftware under the assumption of isothermal 
plasma (ITsuneta et al. Ill99ll ; lHara et al. II1992T ). Only the images without saturated pixels 
in the region of oscillating loops are used. 



3. Method for determining the magnetic field strength based on global 

standing slow mode oscillations 



Applying the MHD wave theory for a straight m agnetic cylindrical model to the coronal 
loop (jEdwin fe Roberts Ill983l ; iRoberts et al. Ill984f ). the oscillation period of a standing slow 



-4- 



mode wave in its fundamental mode is given by, 



„ 2L (I in - 1/2 

P = — , Q = -7 + 



' 4 I 9 1 2 

where L is the loop length, c s = ( r yk B T / [imp) 1 / 2 = 1.52 x 10 4 T 1 / 2 cm s _1 is the sound speed, 
and va = B '/ '(/iop) 1 ^ 2 = 2.18 x 10 n Sne cm s _1 is the Alfven speed. T is the plasma 
temperature (K), B the magnetic field strength (G) in the loop, and n e the electron density 
(cm -3 ). Finally, 7 is taken as 5/3, k B the Boltzmann constant, m p the mass of a proton 
and \x the mean molecular weight, taken here to be 0.6 assuming a standard coronal He 
abundance. 



A recent work of iRoberts I (120061 ) shows that the stratification will modify the period of 
slow modes by, 

21 1 

P= TToi L c « 2vrA, (2) 

ct (1 + L 2 /L 2 ) 1/2 

where A = k B T / fim p g w 50 Mm (T/1.0 MK) is the pressure scale height. For hot loops 
of T=6 MK, it is shown that L c = 1900 Mm, far greater than the observed typical loop 
length. Hence gravitational effects on the period, P, of slow modes in hot SUMER loops are 
negligible. 



From equation (pp) we have 



P= ^ = *±( 1 + l fi f, ,= (1)4. (3 ) 

This equation shows that when the plasma (3 is small, the oscillation period is essentially 
independent of the magnetic field, which is normally the case in the corona. Then the 
derivation of the magnetic field will be very sensitive to uncertainties in period, loop length 
and sound speed measurements. However, in the case of hot loops, j3 is not so small. If we 
assume a coronal loop with an electron density of 10 9 cm -3 and a magnetic field of 10 G, then 
/?=0.06 for cool (1 MK) loops while f3= 0.35 for hot (6 MK) loops. Thus it is reasonable to 
apply this method to hot loops. From equation (pQ) or ([3j), the magnetic field can be derived 
by 

\ 1/2 / r>2 1 \ -1/2 

n 9 \ 1 ( P z 1 x 



where G\ = 4.8 x 10 3 and C2 = 2.3 x 10 4 are constants. B, rig, Tq, P, and L are in units of 
G, 10 9 cm -3 , MK, s and km, respectively. When rig, Tq, P, and L are measured from the 
coordinated SUMER and SXT observations, B can be determined using equation (j4l) . 
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Results 



4.1. SUMER Observations of Hot Loop Oscillations 



We first analyze the Doppler shift oscillations s een by SUMER with the method used in 



Wane et al. 


2002. 


2003a 


b) 



our previous studies (e.g. 
series of intensity and Doppler shift in the Fe XIX line. We find that in all cases the regions 
that exhibited prominent Doppler shift oscillations are located at the loop top. In contrast, 
when the SUMER slit intersected with the loop legs, no noticeable evidence for Doppler 
shift oscillations is found. This is shown in Fig. [7J The southern loop system is the same as 
that shown in Fig. [TJ but closer to the limb. Two intensity peaks occurred simultaneously at 
the legs where the slit intersected with the loop. The time profile of Doppler shifts reveals 
flows in the loop, and only very weak oscillations. Figure [2] shows another example. In this 
case, the slit intersected one leg of a large loop in the north of the analyzed oscillating loop 
(Fig. [2]a). Similarly, the evolution of Doppler shifts shows flows (up to 40 km s -1 ) in the 
loop, but without oscillations. The dominance of the oscillation at the loop top over the 
leg regions is consistent with the features of a standing slow wave in its fundamental mode, 
which has an anti-node at the loop apex and nodes at the footpoints of the loop in velocity 
profiles. 



2005 



Tsiklauri et al. 



Tarovan et al. 120051 ) and the second harmonic mode (e.g. iNakariakov et al. 



1 mode (e.g. 


Selwa et al. 


Nakariakov et al. 


2004; 



20041 ) of standing slow waves can be generated in a coronal loo p. But so far 



no evi dence for the second harmonic has been found in SUMER observations (IWang et al. 
2003bl ). The period of the second harmonic is half the period of the fundamental mode 



and its velocity oscillation has a node at the apex and two anti-nodes at the legs of the 
loop. Therefore, strong (weak) anti-phase Doppler shift oscillations would be expected along 
the two legs (the top). We note that in three cases (Loop-1 shown in Fig. [TJ Loop-3 in 
Fig. [21 and Loop-5 in Fig. 0J the loop top brightenings seem to be associated with Doppler 
shift oscillations in anti-phase, suggesting the possibility of the second harmonics. However, 
we find that the measured oscillation periods are consistent with those estimated for the 
fundamental mode, which are approximately twice the sound wave transit time along the 
loop given the sound speed listed in Table [3] (see Sect. I4.2p . Therefore, this signature is 
probably related to contamination by oscillations from nearby (or super-imposed) loops. For 
example, an arcade system consists of several loops with similar lengths, but their oscillations 
are excited in anti-phase. This explanation is well supported by the case of Loop-1 (the 
northern loop system shown in Fig. [IJa). We find that the southern part of the Doppler shift 
oscillations is more regular and coherent along the slit, while the northern part is complicated 
(Fig. [He). The evolution of intensity along the slit reveals two close brightenings (Fig. [T]&), 
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which correspond to different parts of the oscillations, respectively. In addition, Figure [7] 
also supports that Loop-1 is a loop (or arcade) system. In the case of Loop-3, the two parts 
of the oscillations show different decay times. The southern part is visible for about five 
periods, while the northern part is only visible for about two periods (Fig. [2]c). In the case of 
Loop-5, the two neighboring oscillations are not exactly in anti-phase. The oscillation period 
of the southern one is smaller than the northern one (Fig. H]c). Therefore, events studied in 
this paper are all considered to exhibit fundamental modes. 

We measure physical parameters of the oscillations averaged along the slit in a region 
where the Doppler shift oscillations are strongest (the cuts marked with two parallel lines 
in panel (c) of each figure). For Loop-1, Loop-3 and Loop-5, we choose the regions of more 
regular and coherent oscillations. The function 

V(t) = V + V m em(wt + 0)e~ A ', (5) 

was then fit to the oscillation (Figs. (HEK), where V is the background Doppler shift, V m 
is the Doppler shift amplitude and u, 0, and A are the frequency, phase, and decay rate 
of the oscillations, respectively. We derive the maximum displacement amplitude by A = 
V m /(uj 2 + A 2 ) 1 / 2 . The obtained parameters of the time series are listed in Tabled] We find 
that the oscillation periods in the range of 8—18 min and the decay times of 7—20 min are 



all consistent with the results by IWang et al.l (j2003bl ) in a previous study of 54 oscillation 



events. The oscillation events in Loop-3 and Loop-5 are unusual because they are visible for 
more than 5 period s (Fig. [2]g? andHH). Loop-3 exhibits the slowest decay we have ever found 



(IWang et al. 



2003a] ). with a ratio of the decay time to the period being about 2.4. It was 



associated with a C3.1 flare. 



4.2. SXT Observations of Oscillating Loops 

We measure geometrical parameters of the oscillating loops seen in SXT images, es- 
pecially the loop length, which is a critical parameter to identify the oscillation mode and 
evaluate the magnetic field i n coronal loops as will be shown in Sect. 14.31 We use a method 



similar to that employed by lAschwanden et al.l (120021 ). By assuming a circular loop shape 



this method optimizes two free parameters (ho and 9) to obtain a good match with the ob- 
served loop, where ho is the height of the circular loop center above the solar surface and 9 
is the inclination angle of the loop plane to the vertical. This is shown in Fig. [T]a. We notice 
that two loops, named Loop-2a and Loop-2b, can be discerned in the southern loop system. 
For Loop-2b, an elliptical shape must be assumed in order to get the best fit, in which an 
additional free parameter of ellipticity, e, is included. For Loops 4—7 which were just located 
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above the limb, because of an uncer tainty in determinin g; the positions of their footpoints, we 



apply another method described by lWang et al.l (l2003bl ) in Appendix A. Based on a circular 



model, the geometric parameters of the loop can be derived from measurements of the limb 
footpoint separation and the apex position of the loop by optimizing the free parameter ho- 
Note that with this method the analyzed loop is assumed to have the midpoint between its 
footpoints exactly above the limb. The obtained loop length, the inclination angle and the 
azimuth angle of the footpoint baseline are listed in Table [21 

We compare the evolution of intensity measured at the loop apex by SUMER and SXT. 
Figures [Ul6]e show that the time profiles of soft X-ray flux in the SXT AlMgMn filter and the 
SUMER Fe XIX intensity coincide well, indicating that the Fe XIX obser vations are a good - 



reflect ion of the coronal soft X-ray emission. In our previous studies (e.g. IWang et al.l 12002 



2003al lbL we assumed that hot oscillating loops have a high temperature, 6—8 MK, because 



the Doppler shift oscillations are often observed in the Fe XIX and Fe XXI lines. In this study, 
we directly measure the temperatures of hot oscillating loops by analyzing coordinated SXT 
observations taken through two filters. In all but one event, SXT observations are obtained 
around the times of o scillations (see Fi gs. [LBi and (Ul6]e). With the filter ratio method 



(ITsuneta et al. I Il99ll ; lHara et al. 1 119921 ) , we derive the average electron temperature and 
emission measure for the whole loop. Figures [Ul6]f and [T]l6]g show the evolution of the average 
loop temperature and electron density. The electron density is calculated from the emission 
measure by assuming the plasma depth along the line of sight to be 10 Mm (typical SXT loop 
width) and a filling factor of unity. Some loops, for example Loop-5 (Fig. Hj), show significant 
intensity variations along the loop. This may be because this is a loop system, which is wider 
than 10 Mm and the density could be over-estimated. We find that the temperature during 
the oscillations is in the range of 5—8 MK with a small variation. Most studied loops (except 
for Loop-2 and Loop-4) exhibit a gradual cooling from 7—8 MK and the electron density 
appears to vary slower than the temperature (e.g., Loop-5 in Figs. 0]f and \Qg) . The time- 
averaged temperature and electron density of the loops during the oscillation are listed in 
Table [31 The errors are the standard deviations from the average. Note that since SXT 
missed observations of the oscillation during the main phase for event 7 (Fig. [6]f and ®g) and 
was observing during the decay phase, we take the first data point recorded in the later decay 
phase as the average value. This seems reasonable because the temperature and electron 
density both show a decreasing trend. 
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4.3. Determination of Coronal Magnetic Field 

4-3.1. Magnetic field measurement 

We have measured the oscillation period, loop length, temperature and electron density 
for seven oscillating loops. With these parameters we can calculate the sound speed (c s = 
l^ksT/iimp] 1 ^ 2 ), the tube speed for the fundamental mode of a standing slow mode wave 
(c t = 2L/P), the Alfven speed (va = c s c t /[c 2 s — c 2 } 1 ^ 2 ), plasma j3 (= [2/7] [c 2 /v A ]), and finally 
the mean coronal magnetic field in the loops (eq.jlj). These parameters are listed in Table [31 
and the average values and ranges of these parameters are given in Table [5j If we estimate 
the period of the fundamental mode as ~ 2L/c s , then its ratio to the measured period is 
in the range of 0.75—0.94 for all seven cases, which strongly supports the scenario that the 
loop oscillations we have detected are the fundamental mode. 

We find Alfven speeds in the range 442—1123 km s -1 w ith a mean of 830±223 km s _1 , 



consistent with the result by lOfman fc Aschwanden I (120021 ) who deduced Alfven speeds in 



the range 384—1640 km s -1 with a mean of 986±385 km s _1 for 11 kink-mode oscillating 
loops observed by TRACE. In the case of Loop-2, we have noticed that the loop system 
consists of two loops, both of which intersected the slit at the region where the Doppler shift 
oscillations were detected. For one of these two loops, however, we find that q = 2L/P = 501 
km s _1 , larger than c s . This is inconsistent with ct < c s implied by equation ([1]). This could 
be caused by the uncertainty of the loop length measurement (e.g., neither the circular nor 
elliptical shape is a good approximation for this loop). It is unlikely to be an error in the 
loop temperature because the Fe XIX line, at which loop oscillations are detected, is formed 
within a narrow temperature range. 

We find the plasma (3 in the range 0.15—0.91 with a mean of 0.33±0.26. When excluding 
Loop-4 with (3 = 0.91, (3 is in the range 0.15—0.33 with a mean of 0.24±0.08, which are 
expected for hot (>6 MK) coronal loops. 

Finally, we obtain the magnetic field in the oscillating loops in the range 21—61 G 
with a mean of 34±14 G. These loops have an average height of 32±10 Mm in the range 
24—52 Mm, which are ca lculated by using h = (h + r)cos6 (see Table [2]). In comparison, 



Aschwanden et al.l (120021 ) determined the magnetic field to be in the range 3—90 G from 



TRACE obs ervations of kink-mode o scillat ing loops. With two cases of kink oscillations seen 



by TRACE. iNakariakoy &: Ofman I (120011 ) inferred a narrowed range of 4—30 G. Recently, 
Verwichte et al. I (120041 ) analyzed transverse oscillations of 9 coronal loops of a post-flare 
loop arcade observed by TRACE and interpreted them as standing fast kink oscillations in 
both the fundamental and second harmonics. For those loops, they obtained the magnetic 
field in the range of 9—46 G. Our measurements are basically consistent with these results. 
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On the other hand, it is likely that our measurements over-estimate the magnetic field in 
some events due to an over-estimate of the loop density. Measurements of TRACE loops 
with the temperature of 1—2 MK s howed the electron density to be about (1— 2)xl0 9 cm -3 
(e.g. lAschwanden et al.lll999l . l2000l ). Our measurements give an average electron density of 
around eight times greater. The over-estimate of the loop density could be due to strong 
background emissions in the low corona, which will be discussed in Sect. 14.3.31 



4-3.2. Error analysis 

Now we analyze the uncertainties in the estimated magnetic field. From equation (JJ|) 
the following formula can be derived with the error propagation law, 



and 




Uncertainties in measuring the loop density, temperature, length, and oscillation periods 
all contribute to the uncertainty in measurements of B. Equations (17bl) - (17dj) show that 
contribution by uncertainties in density measurements is the least significant, as f n is a 
constant small number, fx, fi and fp all increase as plasma (3 decreases. Therefore, when (3 
is small, uncertainties in the magnetic field measurement will be very sensitive to the errors 
in the determination of the temperature, loop length and oscillation period. The relations 
in equations (I7b]) - (l7d|) also imply that for relatively high (3 loops, i.e., coronal loops with 
high temperature, high density, and weak magnetic field, there is less relative uncertainty in 
the estimated magnetic field, because of a stronger coupling between the perturbations of 
magnetic pressure and thermal pressure. 

The errors of oscillation period measurement listed in Table [1] are given by the least- 
square fit. The errors of temperature and electron density measurements listed in Table [3] 
are the standard deviations of the time average during the oscillation. Estimates of the 



-10- 



uncertainty for the loop length measurement is difficult, because it involves uncertainties in 
resolving fine loop structures, the loop shape, and locating loop footpoints. For convenience, 
we assume the relative error of loop length measurement io^jV) as 5%, which is comparable 
to those of temperature and electron density measurements (o~x/T and a n /n). By applying 
equation ([6]) the errors of the magnetic field B are calculated and listed in Table [31 For the 
7 analyzed loops, we obtain the relative error of the magnetic field B to be 40% on average 
with f T w 2.4 and = fp ~ 5.9. Given the same set of fx, fh and fp values, the relative 
errors of n e , T, L and P have to be within 10%, 2%, 1% and 1%, respectively, in order to 
improve the accuracy of B measurement to 10%. This suggests that the improvement of 
measurement accuracies in T, L and P is essential. In addition, we estimate the uncertainty 
in L based on the inequality, c t < c s . Since all 7 cases show 2L b s /P = 2(L + o~l)/ P < c s if 
only considering the uncertainty in L, we can obtain a L /L < (1 + 7/3/2) 1 / 2 - 1. Given the 
upper limit of (3 = 0.33 ± 0.26, yielding o"l/L < 22%. This implies that the uncertainty in 
L in our measurements is not larger than 22%, e.g. caused by a departure of loops from the 
circular shape we have assumped. 

We note that the mean density and temperature are used to determine B, while T and n 
both evolve with time because of loop heating and cooling. For more accurate determination 
of B, not only measurements of T and n should be conducted with improved accuracy, but 
also their temporal variations should be delineated from the uncertainty estimates. Theo- 
retically, the evolving T and n may lead to a time dependent loop oscillation period, given 
that the coronal magnetic field and loop length are nearly constant (or their evolution time 
scales are much longer than the cooling time) during the oscillation. Hopefully, high cadence 
observations of time series of Doppler shift oscillations by future missions will allow us to 
examine the time profiles of the oscillation period, P(t), with respect to T(t) and n(t). If 
variations of these three parameters follow equation (j4j), we can minimize uncertainties in 
determining B caused by using the time averaged temperature and electron density. 



4-3.3. Effects of the Background Subtraction 

In addition, we consider the effects of the background plasma emission. Plasma emis- 
sions outside oscillating loops along the line of sight could cause an over-estimate of the 
electron density in the loop (assuming a filling factor of unity). However, to make a back- 
ground subtraction for the whole loop is not an easy matter. An active region loop is rarely 
an isolated loop, and often overlaps with nearby loops along the line-of-sight. This may cause 
a local brightening along the loop (e.g., Loop-5 shown in Fig.HJ). In this case we have no way 
to correct the emissions from other loops. Here we only consider effects of mean background 
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emission in regions adjacent to the oscillating loops in the FOV of the observations (e.g., 
a small box between Loop-1 and Loop-2 in Fig. [T]a). With this mean background emission 
subtracted, we find a slight (~ 3%) increase in temperature in most cases and a decrease in 
electron density for all loops by about 12% (e.g., Figs. [Ijf and d^). We list the time-averaged 
temperatures and electron densities, the derived sound speeds, Alfven speeds, plasma (3 and 
magnetic field in Table H] and their averages and ranges in Table [51 On average, the back- 
ground subtraction causes a decrease in the Alfven speed by about 6%, an increase of the 
plasma (3 by about 15% and a decrease in the magnetic field by about 9%. The effect is 
not significant compared with uncertainties discussed above, specifically uncertainties in the 
loop length, as background emission mainly leads to uncertainties in electron density, which 
has a relatively small influence on B compared to the other parameters. 

We note that the above-mentioned background emission is obtained from a region above 
the oscillating loop. Emission in regions below the loops is much stronger (e.g., the cases 
of Loops 5, 6 and 7). For a comparison, we take the background emission from a small 
box below the loop as the upper limit of the background emission. We repeat the above 
calculations and find that on average the temperature rises by about 17% and the electron 
density reduces by about 40%. Both lead to a reduction of the magnetic field. The combined 
effect leads to an estimate of the magnetic fields with a mean of 22±13 G in the range 
12—51 G. The magnetic fields are reduced by about 35% compared to the values obtained 
before the background subtraction. Therefore, an accurate measurement of coronal loop 
temperature and density is crucial to the determination of magnetic field. 



5. Conclusions and Discussions 

This study presents the first effort to determine mean magnetic field in coronal loops ex- 
hibiting standing slow mode oscillations. We analyzed seven Doppler shift oscillation events 
observed by SUMER in the Fe XIX line that have coordinated SXR imaging observations 
obtained by SXT in at least two filters. This has allowed us to measure the oscillation peri- 
ods of identified loops, and determine the loop geometrical and physical parameters such as 
the loop length, plasma temperature and electron density. Magnetic field is calculated from 
the loop plasma and oscillation parameters. We summarize our results as follows: 

1. It is evident that all oscillation events in this study exhibit the fundamental mode of 
standing slow waves, as pronounced Doppler shift oscillations were located at the loop apex, 
and the oscillation period is approximately twice the sound wave transit time along the loop. 

2. With the filter ratio method, we obtain from the SXT images the temperature and 
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electron density of the oscillating loops. The temperatures are in the range of 5—8 MK, 
consistent with the fact that the Doppler shift oscillations are preferentially observed in the 
hot Fe XIX line. The characteristics of temperature and electron density evolution agree 



with results from numerical simulatio ns with an impulsive heating function (I Warren et al. 
20031 : IPatsourakos fe Klimchuk Iboos h. 



3. Applying MHD wave theory of standing slow mode waves, we obtain for the seven 
oscillating loops Alfven speeds in the range 442—1123 km s _1 with a mean of 830±223 km s -1 
and plasma /3s in the range 0.15—0.91 with a mean of 0.33±0.26. The mean magnetic field 
in the oscillating loops is estimated to be in the range 21—61 G with a mean of 34±14 G. 
Determination of B is sensitive to uncertainties in measurements of temperature, loop length 
and oscillation period. The magnetic field can be derived with higher accuracy in loops with 
a higher (3. Background subtraction in these events could lead to a decrease in the derived 
magnetic field by 9%— 35%. 

Measurements of plasma properties in coronal loops will be greatly improved with the 
upcoming space missions. Solar TErrestrial RElations Observatory (STEREO) will directly 
measure the loop geometry in 3-D, yielding a more accurate loop length measurement. The 
EUV Imaging Spectrometer (EIS) on Solar-B is capable of quick raster scanning to observe 
the Doppler shift oscillation of active region coronal loops both off the limb and on the 
disk. The two EUV bands, 170—210 A and 250—290 A contain many coronal and flare 
lines formed in a wide temperature range, 1—20 MK, which can provide both temperature 
and density diagnostics using the line pair ratio method and a differential emission measure 
(DEM) analysis with several lines. Some studies have indi cated that active region loops may 



have multi-thermal components along the line of sight (e.g. lHara et al. II 19921 : ISchmelz et al. 



200ll ; ISchmelz fe Martens II2006T). so temperatu re measurements through filter and line ratios 
are possibly inadequate (IMartens et al. II2002I ). Instead, we should use the DEM- weighted 
temperature in our analysis. An impulsively generated flare disturbance could excite both 
global fast kink-mode an d slow-mode oscillations in coronal loops as shown by MHD simula- 
tions (ISelwa et al. 1 120061 ) . With joint observations by the X-Ray Telescope (XRT) and EIS 
on Solar-B such events in hot coronal loops could be captured, allowing a direct comparison 
of the estimated magnetic field from both kinds of oscillations. Although the slow mode 
wave has weak connection with the coronal magnetic field compared to the fast kink mode 
wave, an evident advantage is that the occurrence rate of the slow mode o scillations is much 
higher than that of the fast kink mode oscillations (e.g. IWang et al.ll2003bl ). which may make 
the former a more useful diagnostic tool. 



From the theoretical point of view, it is desirable to develop a non-uniform, cylindrical 
magnetic loop model including curvature and stratification effects (e.g. with longitudinal 
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temperature, density and magnetic structures) in order to describe the MHD wave properties 
in a more realistic model, thus can guide or be adapted to the new discoveries likely to emerge 
from new observations, e.g., the varied wave properties along the loop and high harmonics. 
These new observations and theories will lead to significant progress in coronal seismology 
and help us obtain the distribution of magnetic field and other physical parameters along the 
loop, which may serve to solve problems such as coronal heating and solar wind acceleration. 

We would like to thank the referee, Valery Nakariakov, for his valuable comments and 
suggestions. SUMER is financially supported by DLR, CNES, NASA and the ESA PRODEX 
programme (Swiss contribution). SOHO is a project of international co-operation between 
ESA and NASA. Yohkoh is a mission of the Institute of Space and Astronautical Science 
(Japan), with participation from the US and UK. This work was supported by NASA grant 
NNG06GA37G. 
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Fig. 1. — Oscillation events No.l and No. 2 on 19 September 2000. (a) The oscillating soft 
X-ray loops (outlined with diamonds) fitted with a circular model (for Loop-1 and Loop-2a) 
and an elliptical model (for Loop-2b). The SUMER spectrometer slit position is indicated 
as a vertical black box. (b) Line-integrated intensity and (c) Doppler shift time series in the 
Fe XIX line at a fixed slit position, (d) Average time profiles of Doppler shifts along the 
cuts shown in (c). The thick solid curves are the best fit to a damped sine function, (e) 
Average time profiles of line-integrated intensity along the cuts shown in (b). The overlaid 
data points (crosses and diamonds) are the soft X-ray intensity (in arbitrary units) of SXT 
AlMgMn filter averaged in the white boxes shown in (a), (f) Temperature and (g) electron 
density light curves (solid lines) calculated for the whole loop (outlined with the dotted 
curves in (a)) with the filter ratio method. The dotted curves represent the case after the 
subtraction of the background emission, which is taken as the average in a black box, marked 
BK. 
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Fig. 2. — Oscillation event No. 3 on 18 September 2000. 
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Fig. 4. — Oscillation event No. 5 on 17 September 2000, 16:40. 
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Fig. 5.— Oscillation event No.6 on 17 September 2000, 18:30. 



- 21 - 




200 SXT AIMg 
18:07:39 



180 




1 






ij 




- 







900 920 940 960 980 1000 1020 1040 
Solar X (arcsec) 



17:3018:0018:30 17:3018:0018:30 
Time on 2001 Mar 9 (UT) 



-30 I 

(km s~ ) 



CO 



20 



-20 



Q- -40 
o 

Q 





















: (d) V 







17:30 



18:00 



18:30 




^ 6.4 
— ' 6.2 

CD 

-5 6.0 
o 

<d 

^ 5.8 

CD 

^ 5.6 

^ 17:00 
£ 6.0 f 



2 5.5 



SXT A112/AlMg ^ff) ^ 

I 8:00 



17:30 18:00 18:30 

Time on 2001 Mar 9 (UT) 



5.0 
4.5 



° 4 

CD 

^ 17:00 




SXT A112/AlMg (g) 

18:00 

Time on 2001 Mar 9 (UT) 



Fig. 6.— Oscillation event No. 7 on 9 March 2001. 
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Table 1. Time series analysis of Doppler shift oscillation? 



Loop 


to 


v m 


P 







A 




No. 




(km s _1 ) 


(min) 


(rad) 


(min) 


(Mm) 




1 


03:53:53 19-Sep-00 


21 


12.9±0.2 


2.73 


18.1±3.1 


2.64 


4 


2 


04:11:54 19-Sep-00 


36 


18.3±0.7 


-1.47 


19.9±4.2 


6.28 


2 


3 


03:59:25 18-Sep-00 


14 


8.1±0.1 


2.56 


19.8±5.0 


1.05 


5 


4 


04:17:25 17-Sep-00 


17 


8.5±0.2 


-3.00 


7.4±1.7 


1.33 


4 


5 


16:40:09 17-Sep-00 


8 


9.3±0.3 


0.12 


19.1±9.2 


0.68 


5 


6 


18:33:39 17-Sep-00 


30 


8.1±0.3 


-0.66 


7.2±1.8 


2.30 


3 


7 


17:39:37 09-Mar-01 


61 


13.3±0.7 


-3.04 


8.1±2.3 


7.44 


2 



a to is the start time of the modelled time series, V m the Doppler velocity amplitude 
derived by the best fit of a damped sine-function to the oscillations, P the oscillation 
period, the phase of the oscillation, Td the decay time (rj = 1/A), A the displacement 
amplitude (defined as A = V m /(uo 2 + A 2 ) 1 / 2 ), N P the number of periods over which an 
oscillation was detected. 
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Table 2. Geometrical parameters of oscillating loops derived from Yohkoh/SXT images? 



Loop 


lo — Iq 


b - b Q 


a 


e 


i> 


ho 


r 


L 


No. 


(deg) 


(deg) 


(deg) 


(deg) 


(deg) 


(Mm) 


(Mm) 


(Mm) 


1 


-80.2 


10.2 


16.8 


17.6 


20.9 


-22.3 


57.8 


134 


2a 


-68.2 


6.3 


-17.0 


-42.1 


25.9 


20.5 


50.0 


199 


2b b 


-68.3 


6.3 


-18.5 


-44.0 


26.7 


35.2 


70.3/46.4 


275 


3 


-81.4 


8.3 


-26.7 


-24.0 


27.0 


11.4 


21.6 


92 


4 


-90 


10.3 


48.8 


2.9 


48.8 


0.24 


23.5 


74 


5 


-90 


10.5 


18.0 


7.0 


18.0 


-16.9 


42.3 


98 


6 


-90 


8.0 


33.0 


-30.2 


33.0 


-1.9 


27.2 


82 


7 


90 


16.5 


57.9 


35.5 


57.9 


9.0 


37.0 


135 



% — Iq and 6 — b & are the heliographic longitude and latitude relative to Sun 
center for the midpoint of the loop footpoint baseline, a is the azimuth angle of the 
loop baseline to the east-west direction. 6 is the inclination angle of the loop plane 
to the vertical, ip is the angle between the loop baseline and line-of-sight. h is the 
height of the circular loop center in the loop plane, r is the radius of the circular 
loop. L is the loop length. 

b To get the best fit, an elliptical shape is assumed for this loop. The two values 
in r are the semi-major and semi- minor axis lengths, where the semi- major axis is 
parallel to the solar surface. 
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Table 3. Temperature, electron density, sound speed, tube speed, Alfven speed, plasma (3 
and the mean magnetic field inside oscillating loops. 



Loop 


T 


n e 


c s 


Ct 




P 


B c 


No. 


(MK) 


(10 9 cm" 3 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 




(G) 


1 


6.6±0.1 


5.2±0.1 


390 


346 


749 


0.33 


25±6 


2 a 


6.5±0.3 


4.3±0.5 


388 


362 


1025 


0.17 


31±16 


3 


7.0±0.3 


14.1±0.8 


402 


379 


1123 


0.15 


61±30 


4 


6.4±0.9 


10.4±2.0 


385 


290 


442 


0.91 


21±4 


5 


6.3±0.9 


11.1±2.3 


382 


351 


900 


0.22 


43±25 


6 


6.3±0.5 


7.8±0.9 


382 


337 


723 


0.33 


29±10 


7 


5.9±0.3 b 


5.6±0.5 b 


369 


338 


845 


0.23 


29±14 



a In calculations of q, v a, (3 and B, the loop length was taken from Loop-2a in 
Table El 

b In this case T and n e are taken as the maximum values because SXT missed 
observations of the event during the main phase. 

°In calculations of the uncertainty in B, the uncertainty in the loop length is 
taken as 5%L. 
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Table 4. Same as Table [3] but for the cases with the background subtraction. 



Loop 


T 


n e 


C s 








B a 


No. 


(MK) 


(10 9 cm- 3 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 




(G) 


1 


6.8±0.1 


4.6±0.1 


396 


346 


711 


0.37 


22±5 


2 


7.0±0.5 


3.4±0.6 


402 


362 


837 


0.28 


22±9 


3 


7.0±0.4 


12.4±0.9 


402 


379 


1123 


0.15 


57±29 


4 


6.5±1.0 


9.3±2.3 


388 


290 


438 


0.94 


19±4 


5 


6.2±1.1 


10.0±2.5 


378 


351 


943 


0.19 


43±31 


6 


6.7±0.6 


6.8±1.0 


393 


337 


656 


0.43 


25±7 


7 


6.1±0.3 


5.2±0.6 


375 


338 


781 


0.28 


26±10 



a In calculations of the uncertainty in B, the uncertainty in the loop length is 
taken as 5%L. 



Table 5. Average and range of physical parameters of 7 oscillating loops for the cases 
before and after subtractions of the background emission (hk)- 



Parameter 


Before 1^ 


subtraction 


After Jfcfc subtraction 


Average 


Range 


Average Range 


Oscillation period P (min) 


11.2±3.8 


8.1-18.3 


same 


Loop length L (Mm) 


116±44 


74-199 


same 


Temperature T (MK) 


6.4±0.3 


5.9±7.0 


6.6±0.4 6.1-7.0 


Electron density n e (10 9 cm -3 ) 


8.4±3.6 


4.3±14.1 


7.4±3.3 3.4-12.4 


Sound speed c s (km s _1 ) 


385±10 


369-402 


391±11 375-402 


Alfven speed Va (km s~ x ) 


830±223 


442-1123 


784±218 438-1123 


Plasma (3 


0.33±0.26 


0.15-0.91 


0.38±0.27 0.15-0.94 


Magnetic field B (G) 


34±14 


21-61 


31±14 19-57 



